Levy Flights in Inhomogeneous Media 
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We investigate the impact of external periodic potentials on superdiffusive random walks known as Levy 
flights and show that even strongly superdiffusive transport is substantially affected by the external field. Unlike 
ordinary random walks, Levy flights are surprisingly sensitive to the shape of the potential while their asymptotic 
behavior ceases to depend on the Levy index fi. Our analysis is based on a novel generalization of the Fokker- 
Planck equation suitable for systems in thermal equilibrium. Thus, the results presented are applicable to the 
large class of situations in which superdiffusion is caused by topological complexity, such as diffusion on folded 
polymers and scale-free networks. 

PACS numbers: 05.40.-a, 61.41.4e, 02.50.-r, 45.10.Hj 



43 

o 



I 

C 

O 
o 



(N 
O 

c3 



-a 
c 

o 
o 



X 



Diffusion processes are ubiquitous in nature. A freely dif- 
fusive particle is characterized by a mean square displacement 
which increases linearly in time, (X 2 (t)) oc t. However, 
a variety of interesting physical systems violate this tempo- 
ral behavior. For example, the position X(t) of a superdif- 
fusive particle heuristically evolves as X(t) ~ t 1 ^ with 
< /i < 2. Superdiffusion has been observed in a number of 
systems ranging from early discoveries in intermittent chaotic 
systems [1], fluid particles in fully developed turbulence |2], 
to millennial climate changes |3] bacterial motion [4] and hu- 
man eye movements . 

Among the most successful theoretical concepts which 
have been applied to superdiffusive phenomena is a class of 
random walks known as Levy flights 00]. In constrast to or- 
dinary random walks, the displacements Ax of a Levy flight 
lack a well defined variance, due to a heavy tail in the sin- 
gle step probability density. Levy flight models have paved 
the way towards a description of superdiffusive phenomena in 
terms of fractional Fokker-Planck equations (FFPE) 1 8 ] . 

Since many of the aforementioned systems evolve in inho- 
mogeneous environments, it is crucial to understand the in- 
fluence of external potentials on the dynamics. While in or- 
dinary diffusive systems an external force is easily incorpo- 
rated into the dynamics by a drift term in the corresponding 
Fokker-Planck equation (FPE) 1 9] the matter is more subtle in 
superdiffusive systems due to the non-local properties of the 
fractional operators involved. Depending on the underlying 
physical model, different types of FFPEs are appropriate llflll . 
therefore the ad hoc introduction of fractional operators may 
lead to severe problems. 

In cases where the external inhomogeneity can be repre- 
sented by an additive force, considerable progres s has been 
made in a generalized Langevin approach II ll ll 211 which led 
to an FFPE in which deterministic and stochastic motion seg- 
regate into independent components. This approach, however, 
is only suited for systems in which such a segregation can be 
justified on physical grounds. Furthermore, it is only valid for 
systems far from thermodynamic equilibrium and fails in the 
large class of systems which obey Gibbs-Boltzmann thermo- 
dynamics, e.g. where superdiffusion is caused by the complex 
topology on which the process evolves, such as diffusion on 
folded polymers 1 1 311 . 



In this Letter we investigate the impact of external poten- 
tials on this class of systems. Based on the paradigmatic case 
of a randomly hopping particle on a folded copolymer, we 
report a number of bizarre phenomena which emerge when 
Levy flights evolve in periodic potentials and show that ex- 
ternal potentials have a profound effect on the superdiffusive 
transport. This is in sharp contrast to generalized Langevin 
dynamics, which displays trivial asymptotic behavior, a pos- 
sible reason why Levy flights in periodic potentials have at- 
tracted little attention in the past. We demonstrate that even 
strongly superdiffusive Levy flights are highly susceptible to 
periodic potentials. At low temperatures, they exhibit a signif- 
icant dependence on the overall shape of the potential. Coun- 
terintuitively, the asymptotic behavior does not depend on the 
Levy flight index \i, yet differs for various types of potentials 
(except in the ordinary diffusion limit). A perturbation analy- 
sis reveals a universal behavior for high temperatures. Finally 
we show that in finite systems the effect of the potential on 
the generalized diffusion coefficient is least pronounced for 
intermediate values of /i. This is consistent with the obser- 
vation that Levy flights with (j, m 1 are particularly efficient 
in search processes I5l ll4ll . The results we present are a first 
step towards an understanding of superdiffusive dynamics on 
topologically complex structures exposed to external inhomo- 
geneities. 

Let us begin with the dynamics of a particle performing an 
unbiased random walk in a homogeneous environment in con- 
tinuous time. The probability p(x, t) of finding the particle at 
a position x, given that it was initially at the origin is governed 
by the master equation ]9j], 

d t p(x,t) = dy[w(x\y)p(y,t)-w(y\x)p(x,t)], (1) 



in which w(x\y) is the probability rate of initiating a jump 
y — ► x. If this probability has a typical variance in distance 
one may expand the rhs of Eq. ([0 in moments of w{x\y) 
yielding the FPE for a freely diffusive particle, dtp = A p. 
However, when the rate asymptotically follows an inverse 
power law of distance, i.ew(x\y) ~ \x — with /i < 2 

the variance of jump lengths diverges and the particle per- 
forms a superdiffusive walk known as a Levy flight \t(£. In- 
serting this rate into Q the rhs defines the integral operator 
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A^ 2 p(x,i) = J dy 
and Eq. Q may be rewritten as 



\p(y,t) -p{x,t)} 
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d t p = DA^ 2 
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The parameter D is the generalized diffusion coefficient. Up 
to a constant factor, the operator A M ^ 2 is frequently refered 
to as the fractional Laplacian, because it represents a multi- 
plication by — \k\^ in Fourier space |8]. This simple spectral 
property of the operator is the reason why fractional evolution 
equations are frequently introduced in the Fourier domain. 
However, the position representation, Eq. along with 
Eq. 0, provides a more intuitive picture of the dynamics and 
emphasizes the essential fact that A^/ 2 is defined by a non- 
local integral kernel which decreases algebraically with dis- 
tance. Eq. Q is solved by p(x, t) = (Dt) 1 /^ {x/(Dt) x /^) 
where L^(z) = (2iv)~ 1 J dk exp(ikz — |&|^) is the symmet- 
ric Levy stable law of index fi. The argument x/t 1 ^ in re- 
flects the superdiffusive behavior of the process. When /j, = 2 
ordinary diffusion is recovered. In ordinary diffusion, a po- 
tential V is canonically introduced by a drift term (3W'p in 
the FPE. Thus, it may seem reasonable to formally allow for 
an external potential in a superdiffusive system by 

d t p = (3\7V' P + A» ,2 p. (4) 

This type of FFPE has been studied extensively in the past 
Il2ll . It describes deterministic motion in a gradient field 
F = —[3V* subjected to Levy stable white noise r)(t), i.e. 
X = —j3V' + rj(t). However, this approach introduces severe 
restrictions. Systems evolving according to do not obey 
ordinary Gibbs-Boltzmann thermodynamics. The stationary 
state p s , if it exists, is generally not p s oc exp(— (3V) and de- 
pends on the tail parameter fi. Detailed balance is violated, 
and only in the diffusion limit (/x = 2) can the parameter j3 be 
interpreted as an intensive inverse temperature. 

The asymptotics of Eq. in periodic potentials is triv- 
ial. Rescaling the original coordinates x, t — > z = x/j, t = 
t/jV with 7> 1 yields a form invariant FFPE in a new poten- 
tial V(z) — 7 M ~ 2 V(j z). The factor 7^~ 2 <c 1 implies that 
on large spatiotemporal scales any bounded potential is in- 
significant to the dynamics, a possible explanation why Levy 
flights in periodic potentials have attracted little attention in 
the past. 

However, the generalized Langevin description is not ap- 
propriate for a variety of superdiffusive phenomena ficll . as 
a segregation into deterministic and stochastic forces cannot 
be justified by the underlying physics. Consider the system 
depicted in Fig. ^ A particle is loosely attached to a poly- 
mer chain. Thermal activation causes the particle to jump 
between monomers. The heterogeneity of the polymer is ac- 
counted for by the potential V(x) defined on chemical axis x. 
It is reasonable to assume that the rate w(x ± a\x) of mak- 
ing a transition between neighboring sites x and x ± a de- 
creases with increasing potential at the target site and that it 
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Figure 1 : A particle jumps randomly along a folded copolymer which 
consists of two types of periodically arranged monomers. Each 
monomer type has a different potential sketched qualitatively in the 
top left. Due to conformational changes of the chain the particle can 
jump to a site close in Euclidean space yet distant in chemical coor- 
dinate x. 



depends on the potential difference 8V± — V(x ± a) — V(x) 
in units of fc^T = /3 _1 . Both assumptions are accounted for 
by: w(x ± a\x) oc exp[— 05V±/2]. If the polymer is in so- 
lution and subjected to fast conformational changes, regions 
of the chain that are distant along the chemical axis x of the 
polymer may come close in Euclidean space. Long distance 
transitions x — > y may occur with a probability which fol- 
lows an inverse power law of chemical distance \x — y | 
when \x — y\ 3> a. The exponent /x is determined by the 
folding properties of the polymer, e.g. a Gaussian chain im- 
plies \i = 1/2 JlUl . The possibility of initiating distant jumps 
along the chemical coordinate enhances the diffusion process 
considerably, and is believed to play a role in protein dynam- 
ics on DNA strands 1 16]. A jump process of this type can be 
modeled by a master equation Q, in which thermodynamic as 
well as geometric aspects need to be incorporated in the rate 
w(x\y). The above reasoning suggests, 

w(x\y) oc e -Wto- v WV*/\x - y\ 1+ ». (5) 

In Ref. 1 17], a mean field theoretic treatment of the dynamics 
of the polymer chain also yields Eq. 0. Inserting the above 
rate into Eq. Q, we obtain 

d tP = e ^A^ 2 e'^ 2 p-pe^ 2 A^ 2 e' f3V / 2 , (6) 

which is clearly different from the FFPE corresponding 
to generalized Langevin dynamics. Eq. (|6) obeys Gibbs- 
Boltzmann thermodynamics, p s oc exp(— (3V) is the station- 
ary solution, detailed balance is fulfilled, and f3 is a well de- 
fined intensive inverse temperature for all \x £ (0, 2]. Rescal- 
ing coordinates as above yields a potential V(z) = V(~fz) 
lacking the pre-factor 7 2 ~^ which is present in the general- 
ized Langevin scheme. Therefore, the effect of a bounded 
potential will have an effect on all scales. Note that for \i = 2 
Eq. © reduces to the ordinary FPE. When V = the rhs is 
identical to A^/ 2 p. 

Letting ip(x,t) = exp[— (3V(x)/2]p(x, t), Eq. © can be 
recast into a fractional Schrtidinger equation, 

dti> = -n^ (7) 

H = -A^ 2 +U, U= e PV/2 Afi /2 e -0V/2 (8) 

with an anomalous kinetic term — A p / 2 and an effective po- 
tential U which depends on A separation ansatz yields the 
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Figure 2: Bandstructure k„ i? as a function of for four different po- 
tentials. (A) a simple cosine potential V oc cos(a;/A), (B) the square 
wave potential, and potentials given by V oc ±(1 + cos(a;/A)) 7 . 
When 7>1, the latter possess localized high potential barriers (C) 
or localized potential wells (D). A Levy flight (/i = 1 /2) is compared 
to ordinary diffusion (/i = 2). 



associated stationary equation, 

E + A^ 2 - u] ip(x) = 



(9) 



for the spectrum E. Let us consider periodic potentials of 
wavelength 2n\, V{x) = V(x + 2n\n), with bgZ. With- 
out loss of generality we restrict ourselves to potentials with 
vanishing offset and unit variance. A Bloch ansatz tl> q (x) — 
e iqx 0{x) with 0{x) = 6(x + 2n\n) and q e [0, 1/A] inserted 
into |7} in Fourier space yields 



E n,q 



U 



o 



(10) 



In Eq. (II Oi E n . q are the eigenvalue bands labeled by the dis- 
crete band index n and the continuous Bloch phase q. The 
spectrum of the system for vanishing potential is given by 
E® = |n/A — q\^. The Fourier coefficients of the peri- 
odic component of the eigenfunction and the effective poten- 
tial are given by 9 n = l/2ir\ J 2 x dx9(x) exp(— inx/X) and 
U n = l/2%Xf„ , dxll(x) exp(— inx/A), respectively. The 
spectrum E n>q depends implicitly on (3. In the high temper- 
ature limit the eigenvalue bands merge to form a continuous 
spectrum. For non-vanishing j3 gaps between bands emerge. 
The band structure determines the relaxation properties of the 
system. For a comparison of different Levy indices fi, it is 
more appropriate to compare the generalized crystal momen- 
tum defined as k„ i9 = E„{ q . 

Fig.|2]depicts the /3-dependence of the first few bands K n>q 
for different potentials, each one reflecting a frequently en- 
countered physical situation. In each panel ordinary diffusion 
is compared to enhanced diffusion with Levy index \i = 1/2. 
The band structure of ordinary diffusion (n = 2) in the cosine 
potential displays only one significant gap 1 18] contrasting the 




Figure 3: The quantity (? M (A/L) in the high temperature regime as 
a function of Levy index ji and fixed system size (left). The solid 
lines (symbols) depict the results obtained from perturbation theory 
(numerics). Thick gray lines indicate the asymptotic limit of 1/4 (1) 
if (j, < 2 (jj, = 2). Viewed as a function of inverse relative system 
size X/L for a set of values of /i (right) indicates that convergence to 
the limiting values is slowest when /i is small or slightly less than 2. 



superdiffusive case in which the narrowing effect of individ- 
ual bands is substantial (Fig.|2j\). The effect is even more pro- 
nounced in the square wave potential (Fig.|2jj). Band coupling 
leads to far more complex band structures when [i = 1/2. 
The most striking difference occurs in the localized barrier 
(well) potentials, Fig. |2jT(D). In the example of copolymers 
discussed above, these cases describe situations in which the 
polymer consists mainly of a single type of monomer inter- 
spersed with small intervals of another type of monomer at 
a higher or lower potential, respectively. On one hand, the 
band structures are identical when \i = 2, indicating that an 
ordinary diffusion process does not distinguish between bar- 
riers and wells. On the other hand, if \i = 1/2, the band 
structures differ considerably, the shape of the potential has a 
profound impact on the band structure and thus on the dynam- 
ics of the system. In a system containing barriers, a particle 
located at any given position may initiate a distant jump with 
a probability which is decreased by the concentration of en- 
ergetically unfavorable target positions. In the repetitive bar- 
rier potential this is low, so the process is not considerably 
affected. In Fig.|2jD the particle is likely to be trapped in a po- 
tential well. The probability of initiating a distant transition to 
another energetically favorable state is low and the repetitive 
well potential slows down the dispersion dramatically. 

The asymptotic behavior is governed by the lowest band of 
the spectrum, n = 0, and q <C A -1 . For (3 = the system 
is freely superdiffusive. This yields E® = q^. When (3 is 
finite the (^-dependence remains, i.e. E q m D(f3) q^. The 
generalized diffusion coefficient D((3), however, is reduced 
to a value less than unity; the process is slowed down when 
a potential is present. The high temperature regime can be 
investigated by expanding the effective potential U in (|9j in 
powers of (3. Neglecting all terms of order higher than 0(f3 3 ) 
we obtain 



U n 



1 

a7 



Y^Vn-miWT -2\m\^)V„. 



(11) 



Inserted into (II 01 yields in perturbation theory Eo. q w g M [l 
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/3 2 G M (q)]. As expected, the generalized diffusion coefficient 
D((3) = 1 — p 2 G IJi (q) decreases quadratically with increasing 
(3. The factor G M (q) is positive and depends on V and \i. It 
quantifies the effect on the asymptotics, the larger G^(q) the 
stronger the slowing down effect of the potential. We obtain 



^0?) = tX! \ v m\ 2 9^(q^/m) where 



(12) 



m>0 



1 



zY - z» (1 + zY - z^ 



2 . 



Noting that J2 m \ Vm\ 2 — 1 the asymptotic limit is 



HmG^g) = 



< /i < 2. 



(13) 



(14) 



The asymptotics are the same for any type of potential. In ad- 
dition, the rhs of dl4> is independent of fi with a discontinuity 
to a higher value on the margin fi = 2. 

The limit q — > represents an idealized system of infinite 
extent. In a finite finite system of size 2irL, the Bloch phase 
acquires discrete values q = n/L with n G N. The relaxation 
time is defined by inverse of the lowest eigenvalue, obtained 
by d!2t at q = L~ x <A _1 . The result is shown in Fig.|3]for 
the cosine potential. In this case, eq. H2\ implies G p (1/ L) = 
g^X/L). On the left, g^ is depicted as a function of fx for 
a number of system sizes. Surprisingly, the asymptotic limit 
is not attained uniformly on the /^-interval (0,2]. Even for 
very large systems g^ exhibits a minimum at an intermediate 
value /i w 1. Interestingly, as /i — > the factor g^ diverges. 
Although small values of /i are equivalent to heavy tails in 
the transition probability, the potential strongly influences the 
dynamics in that range. We conclude that Levy flights with 
intermediate values of \i are most robust when perturbed by an 
external field. This may explain why Levy flights with (x « 1 
are the most efficient when employed in random search lll4ll . 

Finally, we investigate the effective generalized diffusion 
coefficient D{j3) in the low temperature regime. Since pertur- 
bation theory fails here, we must rely on the numerical diago- 
nalization of dlOi . The result is depicted in figure^] D(j3) is 
identical for all superdiffusive processes in a given potential. 
However, a comparison between potentials reveals a unique 
response of Levy flights to each potential shown in the inset. 
The effect on D{(3) is least pronounced in the potential barrier 
system, intermediate for the cosine, and strongest in the po- 
tential well. In contrast, ordinary diffusion shows a decrease 
in D{(3) which is not only greater compared to all the other 
cases, but is independent of the shape of the potential. For 
small (3 the results are consistent with those obtained from 
perturbation theory as indicated by the dotted (p, < 2) and 
dashed = 2) lines. 

In this Letter we have shown that Levy flights are substan- 
tially affected by external inhomogeneities, a feature generic 
to topologically induced superdiffusion, absent in popular 




Figure 4: The generalized diffusion coefficient D(ff) as a function 
of inverse temperature for a chosen set of Levy indices (indicated by 
the symbols in the lower left) and the set of potentials depicted in 
the upper right and distinguished by color. The data were obtained 
by numerical diagonalization of The dotted (fi < 2) and dashed 
(/1 = 2) curves indicate the result obtained by the perturbation ex- 
pansion for small (5. 



generalized Langevin models, yet crucial for the understand- 
ing of physical applications such as protein search on DNA 
strands and random motion on complex networks. 

D. Brockmann thanks I.M. Sokolov and W. Noyes for in- 
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